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Starting from the chiral Lagrangian for Wilson fermions at nonzero lattice spacing we have ob- 
tained compact expressions for all spectral correlation functions of the Hermitian Wilson Dirac 
operator in the £-domain of QCD with dynamical quarks. We have also obtained the distribution 
of the chiralities over the real eigenvalues of the Wilson Dirac operator for any number of fla- 
vors. All results have been derived for a fixed index of the Dirac operator. An important effect 
of dynamical quarks is that they completely suppress the inverse square root singularity in the 
spectral density of the Hermitian Wilson Dirac operator. The analytical results are given in terms 
of an integral over a diffusion kernel for which the square of the lattice spacing plays the role of 
time. This approach greatly simplifies the expressions which we here reduce to the evaluation of 
two-dimensional integrals. 
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1. Introduction and Motivation 

The low-lying spectrum of the continuum Dirac operator has been described successfully in 
terms of chiral Lagrangians or Random Matrix Theory The reason for this success is that the 

smallest eigenvalues are of order l/V , which is deep inside the £-domain. In this domain the static 
part of the chiral Lagrangian, which is equivalent to Random Matrix Theory, factorizes from the 
partition function, and determines properties of the low-lying Dirac eigenvalues. In this lecture we 
discuss the extension of these results to nonzero lattice spacing. We can distinguish two different 
types of discretization effects: those that do not affect the symmetries of the Dirac operator and 
those that do affect its symmetries. In the first case, discretization effects can be absorbed into a 
redefinition of the low-energy constants. In the second case, the discretization errors take us outside 
the universality class which requires additional terms in the chiral Lagrangian. The two main 
discretization schemes are staggered fermions and Wilson fermions and both break the symmetries 
of the continuum Dirac operator. In this talk we will discuss recent results [||, 0, |5|, ^ for the the 
Wilson Dirac operator including its dependence on the topological index and the number of flavors. 
The confirmation of these results by lattice simulations [^, ||] will be discussed elsewhere in these 
Proceedings. For recent progress on the staggered Dirac operator see [^. 

We will consider discretization effects of the Wilson Dirac operator 

Oiy = ^7M(VM + V*)-^«V*V^. (1.1) 

Because of its 75 Hermiticity it is advantageous to analyze the spectrum of the Hermitian Dirac 
operator [0 [n|, |l2[ [Tl] 



D5 = y5{Dw + m)=Dl. (1.2) 

Since Dyv is nonhermitian, its eigenvalues are complex, but because of the 75 Hermiticity, they 
occur in complex conjugate pairs or are real. Therefore the complex eigenvalues can collide with 
the real axis and turn into a pair of real eigenvalues. This is different from QCD at nonzero chemical 



potential 1 14 ] where the nonzero Dirac eigenvalues are complex and occur in pairs of opposite sign. 
Then both their real and imaginary parts have to vanish for a pair to collide with the real axis. 

In this talk we report on results for discretization effects on spectra of the Wilson Dirac opera- 
tor of QCD with dynamical quarks. A direct evaluation of the generating function is complicated. 



but using a diffusion method Ill5|] we obtain compact expressions for correlation functions of Dirac 



eigenvalues which we here further simplified by exploiting an underlying Pfaffian structure []16|]. 

We start from a chiral Lagrangian, but our results can also be derived from the corresponding 
chiral Random Matrix Theory [||, §]. Recently, the joint eigenvalue distribution of both the Wilson 



Dirac operator [17] and the hermitian Wilson Dirac operator [18| have been obtained from RMT. 



2. Spectrum of the Dirac Operator 

In order to access the spectrum of D5, we introduce an axial mass z 

det(Z)„, + m + 75z) = det(75(Dvv' + m)+z). (2.1) 
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Figure 1: Spectral flow of eigenvalues of D5 versus the quark mass m. 



The low-energy limit of the corresponding partition function is given by a chiral Lagrangian that 
up to low energy constants is uniquely determined by symmetries 20, 21] 

if (m,z;a) = ^mITr([/ + U'^) + ^zITr([/ - U'^) - a^WgTv{U'^ + U'^). (2.2) 

In the microscopic domain, where the combinations mV, zV and a^V are kept fixed in the thermo- 
dynamic limit, the m, z and a-dependence of the partition function resides in its zero momentum 
part that factorizes from the nonzero momentum part 

Zlim,z;a) = I Jf/det^[/^^mVZTr(f/+J/t)+^.yZTr(t/-J/t)-«^VW/,Tr(f/^+f/-^)^ ^2.3) 

' JueU{Nf) 

To this order the chiral Lagrangian also contains —a^W(,[Tr{U + U^)]'^ — a^Wj[Tr{U — U^)]'^ but 
they can be added to the mass term after linearizing the squares at the expense of a Gaussian 
integral and will not be considered below. 

Since small Dirac eigenvalues are very sensitive to the index, v, of the Dirac operator we 
will work at fixed V. This index is determined from the spectral flow of the eigenvalues of D${m) 
[22, 23] (see Fig. 1). At the crossing point, m^, with the real axis we have (for small a, the physical 
part of the lattice Dirac spectrum can be separated from the unphysical part) 

75(Div + m,.)0 =0 => Dw^ = -mc^. (2.4) 

Note that multiple crossings may occur. Therefore the number of real eigenvalues of D„ is at least 
equal to its index. The total number of flow lines with a net flow across the real axis is a topological 
invariant of the Dirac operator. It can be written as a sum over eigenvalues of D5 [22] 

v= E signf^^) =Esign((0,|75|0.». (2.5) 

mk,X\mk)=0 V /m=mj. k 

The 6 -dependence of the QCD partition function is given by 

ZNf{m,z;a,d) =Y^e'^^Z)^^.{m,z;a). (2.6) 

V 

In the continuum limit, because of the anomaly, the factor exp(jv0) can be written as the Jacobian 
of an axial transformation on the quark fields. At nonzero a the natural extension is (See [24].) 

m ^ e'^^'/^fm, a ^ e'^'^'/^fa. (2.7) 

As is the case in the continuum limit, the partition function at fixed v is defined by inverting the 
decomposition (2^), which also can be applied to the chiral Lagrangian 

Z](i^.{m,z;a) = dee'''^ZNf{m,z:,a,d), (2.8) 
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where the 0-dependence of the chiral Lagrangian is induced by the transformation (2.7). 
The generating functions for the Wilson Dirac spectrum are given by [|3|] 

The resolvent (/? = !) and the spectral density are given by 



, p^(A^,m;<3) = —ImG^ (z,m;a) I 



For z <^ Aqcd the z-dependence of the generating function is given by a chiral Lagrangian that 
is uniquely determined by symmetries that are compatible with the convergence of the bosonic 
integrals. In the microscopic domain the generating function reduces to the graded integral [jsj, ^ 



^uAm,z,a) = , } , / ^^sdet^([/)^.'i'"^^Tr(j/-[/t)+,-iy2:Tr(?[/+j/t?)-/\,2yWsTr(j/^+[/-2) 
^ ' Zl^^{m;a) JueGliNf+i\i) 



(2.9) 

where, (^3 = diag(0, • • • ,0,z,z')- The transformation U — )■ iU is required to get convergent integrals 
in the noncompact sector for Wg > (for which the Aoki phase [25] can occur). 

A direct evaluation of the integrals over the supergroup becomes very laborious for Nf > 1 
but they still can be worked out for one flavor ^ However, for degenerate quark masses the 
expressions can be simplified by starting from the identity 



where a is an {Nf+p\p) graded "Hermitian" matrix. In a diagonal representation of a the partition 
function can be expressed in terms of an integral over shifted quark masses. The result for the 
resolvent is given by 

Gl^^,^,{m,z;a) = — i— ^ / J57(5).(^/^^''^)^^^(^-)'sdet^(5-m)Z^^,^iu({^m2-52};a = 0), 
'^f^ ' ' (2.10) 

where J{S) is the Jacobian of transforming a to the diagonal representation S. 
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2.1 Effect of Light Flavors 

For a = the spectral density of 75 {Dy/ + m) is obtained by a simple variable transformation 
of the spectral density p^{x) of Dyv: pl{x) = x/^/x^ -rrfip'^{\/x^ -m?-)B{\x\ - \m\) (See Fig. | 
(right)). For a^O the result for the spectral density follows from the resolvent ( 2.10| ) (see Fig. § 
(left)). Although, for Nf = 1 and Nf = 1 the integrals in ( 2.10| ) can still be evaluated directly, it 
becomes computationally expensive. However, because of an underlying Pfaffian structure [jl^], 
the integrals can be written as the sum of products of two and one dimensional integrals. 

The three dimensional integral for one flavor resolvent (2.10) can be rewritten as 



1 (z, m; a) = Gi 1 1 (z, z, m; a) + Gi 1 1 (0, z, m; a) \ 



Zl, {z,m;a) 
'ZUO,z,m;a)^ r, (2-11) 



where the first term is the quenched resolvent G^-^{z,z,m;a) with 
The various partition functions that enter in the expression aie given by 



ZQn(z,m;a) 
Zi{z,m;a) 
Zl{z]_,Z2,m;a) 



1 



2_/2w/2 



1 

1 

xZl{{m'+s\)"\{m^ + slYl^), 
Zl,{x,y) = ^-^[yK^+,{y)h{x)+xK^{y)h+,{x)l Zl{x,y) 



(2.12) 

yly+l {y)ly (X) - X/y+l {x)ly (j) 



Likewise the four-dimensional integral ( [2.10| ) for the Nf = 2 resolvent be re-expressed as 

G^n(z,m;a) (2.13) 



G\,Az,z.,m;a) 



idGi\i{zi,z,m;a) 



dzi 



ZX{0,z,m;a) G\,^{0,z,m;a) dZUzuZ,m;a) 
+ 1- 



1=0 



Z2 (m;^) 



ZX{m;a) 



dzi 



7A=0 
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where 

xZl{{m^+s\Yl^,{m' + sl)'l^). (2.14) 
These simplified expressions for the unquenched resolvents have not been presented previously. 

2.2 Density of Real Eigenvalues 

The Wilson Dirac operator has as least v real eigenvalues in the physical sector of the Dirac 
spectrum. Because the eigenvalues of the Wilson Dirac operator occur in complex conjugate pairs, 
additional real eigenvalues may be produced when a pair collides with the real axis. Depending on 
the expectation value of 75 the density of the real eigenvalues may be decomposed as [ p7| ] 

Pr:ai(A^)=p;(A>^) + pAA^). (2.15) 
The distribution of the chiralities over real part of the spectrum of D}/^ 

Pl(^^) = l I 5(A'^ + Ar)sign((fc|75|^») = p^A^) - p;(A«'), (2.16) 

was derived from the chiral Lagrangian {IS) both quenched and with dynamical quarks 

[|,|]. We have that 

Pz(A'^)<Pr:ai(A^) and /^A«'p^^(A^) = ( £ sign((^|75|^))) =v. (2.17) 

The expression for p^(A^) can again be simplified to at most two-dimensional integrals. 

The calculation of the density of real eigenvalues is more complicated. However, starting from 
a random matrix theory for the Wilson Dirac operator, it was possible to derive the joint eigenvalue 



probability density of Dw [ ]17| ] with the total number of real eigenvalues given by [ ]17[ , [27| ] 



3. Conclusions 

We have obtained the lattice spacing dependence of the microscopic spectral density and 
all spectral correlators of the Hermitian Wilson Dirac operator with degenerate dynamical quark 
masses. This makes it possible to extract the low-energy constant from the distribution of the 
smallest eigenvalues. It should satisfy the constraint due to 75-Hermiticity that Ws > (see also 
the recent discussion in [^]). These results can be extended to include the two other low-energy 
constants to this order. Our results also explain the lattice result that the width of the distribution 
of the smallest eigenvalue scales as ajy/W . More results can be obtained using a random matrix 
formulation of the chiral Lagrangian. Then powerful random matrix methods make it possible to 
obtain all spectral correlation functions of D^^^ in the microscopic domain. 
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